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Abstract. An Ising-type classical statistical ensemble can describe the quantum physics of 
i <~| fermions if one chooses a particular law for the time evolution of the probability distribution. It 

Oh accounts for the time evolution of a quantum field theory for Dirac particles in an external 

electromagnetic field. This yields in the non-relativistic one-particle limit the Schrodinger 
equation for a quantum particle in a potential. Interference or tunneling arise from classical 
probabilities. 



1. Introduction 

It is often believed that quantum physics and classical statistical physics are fundamentally 
different on a conceptual level. For both branches of physics probabilities play a central 
role. Classical statistical physics involves only the concept of a positive definite probability 
distribution {p T (t)}, rules for the computation of expectation values of observables and 
correlation functions in terms of {p T (t)}, and a dynamical equation describing the time evolution 
of the probabilities. In quantum physics the complex wave function plays a central role. It 
is often said that quantum physics is the physics of phases in addition to the probabilities. 
Observables and correlations are associated to non-commuting operators. These properties lead 
to characteristic features as interference, entanglement, tunneling, fermions and bosons, that 
are all not familiar in macroscopic classical statistical physics. 

We present here a microscopic classical statistical model that accounts for all characteristic 
features of quantum mechanics. We consider a generalized Ising-type model for discrete Ising 
spins on a lattice. The discrete variables can be associated with simple yes/no alternatives 
or bits. The states r of the classical statistical ensemble are configurations of Ising spins or 
lattice-ordered sequences of bits. The probabilities p T are positive and normalized, p T > 0, 
S T p T = 1. We will propose a particular evolution law for the time dependence of the probability 
distribution {p T (t)} that respects these properties. We will show that for this law the Ising-type 
classical statistical model describes a quantum field theory for Dirac-fermions in an external 
electromagnetic field. Specializing to probability distributions that describe one-particle states 
we recover the Dirac equation in an external electromagnetic field. In the limit of a static electric 
potential and for a non-relativistic particle this yields the Schrodinger equation for the quantum 
wave function in a potential. The characteristic features of quantum mechanics as interference 
in a double slit experiment are therefore accounted for by a suitable evolution law for classical 



probabilities. Quantum mechanics, quantum fermions and quantum field theory can be obtained 
from classical statistics. 

This result has far reaching conceptual, and in the future perhaps also practical, consequences. 
Quantum systems turn out to be a particular class of classical statistical systems, characterized 
by particular properties of the evolution law. They are now embedded in a wider class of 
classical statistical systems with different evolution laws. This opens the possibility to describe 
phenomena as decoherence [1] or syncoherence [2] not only by embedding a quantum subsystem 
into a larger quantum system that includes the environment. Many characteristic features of 
decoherence and syncoherence may find a description in terms of probabilities for the subsystem 
alone, but with an effective evolution law different from the unitary quantum evolution. Precision 
tests of quantum mechanics can be described in a coherent conceptual framework for "zwitter"- 
particles whose properties interpolate continuously between a quantum particle and a classical 
particle [3]. On the other hand, a unified conceptual framework can also help to understand the 
observed quantum features of interference and tunneling in experiments with classical objects, as 
droplets of a liquid [3J. Finally, one may hope that the equivalent classical statistical description 
of quantum systems may permit to device new calculational tools. 

It may be worthwhile to specify clearly what we mean by a classical statistical system. The 
first ingredient is the classical probability distribution {p T (t)}. The second one is the rule for 
the computation of expectation values of observables. Classical observables A have a fixed value 
A T in every state r. These are the only possible outcomes of a measurement. The expectation 
value obeys (A) = Y^ T A T p T . Third, one needs to specify a law for the time evolution of {p T (t)}. 
In our model the Ising spins are the fundamental entities. We do not assume the existence 
of some underlying law that determines the time evolution of a given spin configuration. Our 
concept is "probabilistic realism" , where the description of nature is genuinely probabilistic. The 
basic dynamical law should therefore be formulated as an evolution law for the probabilities. 
Conceptually, this is similar to a formulation of classical mechanics with the Liouville equation as 
basic dynamical equation. (Newton's laws for the motion of planets can, of course, be recovered 
from the Liouville equation for the particular ensemble which describes isolated bodies.) In 
classical mechanics we can also start from Newton's equation and derive the Liouville equation, 
interpreting the probabilities as a "lack of knowledge" of a deterministic system. This possibility 
is not available in our setting since we have no evolution law for individual spin configurations. 
(It is not formally excluded that an evolution law for individual spin configurations exists, 
similar to cellular automata [5]. It will, however, be difficult to find it and we do not need it.) 
Probabilities do not reflect the lack of knowledge of an otherwise more complete system. They 
are simply the basic concept for formulating laws of nature. On a fundamental level we only 
require that the evolution law preserves positivity and normalization of the probabilities and 
that it is causal in the sense that the distribution {p T (t + e)} at time t + e can be computed from 
the distribution at time t. This typically leads to a differential evolution equation. In practice, 
the proposed evolution law will also preserve locality in space. (No signals travel faster than 
light.) 

Our derivation of the Schrodinger equation in terms of an Ising type classical statistical 
ensemble and the description of the outcome of a double slit experiment only involves these 
three basic ingredients of a classical statistical system. In particular, we consider "interval 
observables" , which take the value one if a particle is present in a certain region of space, and 
zero if it is absent. These observables are directly formulated in terms of the discrete classical 
observables for the Ising spins. The "discreteness" of individual particles is therefore a natural 
outcome. The continuous aspects of the wave function will be related to the description of 
probabilities by continuous functions. Particle- wave duality finds a natural explanation in our 
probabilistic setting. We implement the classical statistical rule that only values A T in the 
spectrum of an observable are found in measurements. This distinguishes our approach from 



Bohm's quantum mechanics [6j or other approaches where the wave function is a classical field, 
which has a continuous spectrum and is, in principle, observable. 

The fourth ingredient for a classical statistical system is a rule for the computation of a 
correlation function that describes a sequence of ideal measurements. It should be emphasized 
that such a rule is an independent postulate for any probabilistic description. As every 
experimentalist knows there are many ways to perform a sequence of two measurements - 
"ideal" ones and "less ideal" ones. Correlations describe the outcome of a sequence of two 
measurements. They depend on the precise way in which measurements are performed. "Less 
ideal" measurements will yield correlations which differ from "ideal" measurements. The 
dependence of the correlation function (BA) on the precise way of performing a sequence of 
measurements for A and B is reflected mathematically in the many different possibilities to 
define a product structure between observables. A product structure assigns to two observables 
A and B a third observable C, namely (B, A) — > C = B o A. Translated to experiments, C 
corresponds to the "combined observable" whose expectation value is the correlation function. 
In general, the correlation function (BA) depends on the conditional probabilities w(B p \A T ) 
to find for B the value B p if a measurement of A has found the value A T . These conditional 
probabilities depend on the precise way how the sequence of measurements is performed. 

Since there are many ways to perform sequences of measurements, and many associated 
correlation functions, one wants to define the notion of "ideal measurements". Correlations 
differing from the ones for ideal measurements can then be attributed to "imperfect 
measurements". It is at this point that some thought is needed. For macroscopic classical 
statistical systems an idealized sequence of measurements is generally associated to the classical 
correlation function (BA) = ^2 T p T B T A T . While well adapted to large macroscopic systems, 
as two temperature measurements at different positions in a liquid, it often fails for isolated 
microscopic subsystems as an isolated atom. For an ideal measurement of properties of a 
subsystem the correlation function has to be computable from the information available for 
the description of the subsystem. This requirement often conflicts with the classical correlation 
function [2j. Classical correlation functions often involve properties of the subsystem as well 
as its environment (apparatus). The appropriate correlation functions for ideal measurements 
differ from the classical correlation function in these cases. 

We emphasize that the issue of correlations is not relevant for the description of the time 
evolution of expectation values of observables. Thus our description of interference, tunneling 
or entanglement in terms of the time evolution of a "classical" probability distribution does 
not involve any choice of the correlation function. This choice matters, however, if we want to 
describe the outcome of sequences of ideal measurements. It is important for an understanding 
why Bell's inequalities [HE] are violated for the correlations in quantum systems [2J. 

The present work is based on a series of papers by the author [DJ, [ID], [3], [2J and has 
substantial overlap with ref. [DJ. 

2. Generalized Ising model 

We will consider an Ising-model type system for discrete Ising-spins, s 7 (x) = ±1, with x 
points on a suitable three-dimensional lattice and 7 denoting different "species" of Ising-spins, 
7 = 1 . . . N s . For N s = 4 our system will be equivalent to a quantum field theory for Majorana 
or Weyl spinors, while for N s = 8 we will describe Dirac spinors. The states r are sequences or 
configurations of Ising-spins, r = {s^(x)}. Instead of Ising-spins, we will actually use occupation 
numbers or bits n 7 (x) = (s^(x) + l) /2, such that the states r describe bit sequences of numbers 
n 7 (x) = 0,1, r = {ra 7 (x)|. We will consider a setting with L 3 /8 lattice points (L even) such 
that the number of states is 2( JVsi3 / 8 ). The continuum limit L — > 00 is taken at the end. 

The sequences or configurations of occupation numbers |n 7 (x)} show already a strong 
analogy to the basis states of a quantum theory for an arbitrary number of fermions, formulated 



in the occupation number basis for positions. We will exploit this analogy in order to formulate 
a fundamental "evolution law" for the time evolution of the probability distribution p T (t) , such 
that our system describes a relativistic quantum field theory for fermions. 

2.1. Classical wave function 

In order to define our model we need to specify an evolution law for the time dependence of the 
probability distribution {p T (t)}. A useful concept for this purpose is the classical wave function 
{?r(*)} [3J, defined by 

q T (t) = Sr(i)VP^) . Pr(t) = . S r(t) = ±1- (1) 

This is a real function which is given by the square roots of the probabilities up to signs s T (t). 
In terms of the wave function the expectation values of classical observables obey the "quantum 
rule" 

{A(t)) = (q(t)A q (t))=^2p T (t)A T , (2) 

T 

with A a diagonal operator (A) rp = A T 5 rp . The quantum rule follows directly from the classical 
statistical rule. The signs s T do not appear in eq. (J2J). At this stage the use of a wave function 
seems redundant. Nevertheless, it permits the use of the quantum formalism for a description 
of a classical statistical system. 

Consistent evolution laws have to correspond to rotations of the vector q T (t), 

ir{t) = Yl ^pfo *03p(*o . Yl R ^ o^pfo o = (3) 
p p 

The positivity of p T = is automatic, and the normalization of the distribution remains 
preserved provided it was normalized for some initial time £j n , since the length of a vector 
or X T q^ = T, T p T is preserved by rotations. We will specify a linear evolution for the wave 
function. In other words, we will consider rotation matrices R Tp that are independent of the 
wave function. The classical wave function specifies the probability distribution uniquely. The 
specification of an evolution law for the wave function therefore defines the dynamics of the 
classical statistical system completely. 

Quantum wave functions are usually defined in a complex Hilbert space. Indeed, many 
characteristic quantum features are closely associated to the "physics of phases" . In our setting 
we will define later a complex quantum wave function by introducing a complex structure 
in the real space spanned by {q T }- The 2^ NsL / 8 ) real components of the vector {q T } can 
then be associated to a complex vector with dimension 2^ NsIj3 ^ s ^ 1 . There is no conceptual 
difference between the classical and quantum wave function in our setting. However, more of 
the information contained in a complex wave function ip a = ipR,a + iipi,a wm now directly be 
associated to classical probabilities. In quantum mechanics one usually distinguishes between 
the probabilities p a = \ip a \ 2 an d the phases, whereas in our approach both positive numbers 
ty\ a and V>? a will be associated to classical probabilities p T . Our description of fermions and 
quantum field theory will be based entirely on the real "classical" wave function {q T }, which 
corresponds to the roots of the probabilities according to eq. JT]). 

The only information contained in the classical wave function which goes beyond the set of 
probabilities {p T } concerns the signs {s T }. The choice of {s T } does not affect the expectation 
values of classical observables and may therefore be considered as a type of gauge choice. A linear 
evolution law involves, however, directly these signs and one may wonder if the evolution 
of the classical wave function involves physics beyond the classical statistical setting, thereby 
"smuggling in" quantum mechanics in some hidden way. This is actually not the case. For 



computing q(t + e) from a given wave function at time t according to eq. ([3|), we only have for 
every r the choice between the two values s T (t) = ±1. Only one of them is consistent with the 
normalization of the probabilities. This fixes for all t the sign distribution {s T (i)} in terms of 
some initial distribution" {s T (to)} for t = to. For a differential evolution equation it is easy to 
see [3] that s T (t) is constant as long as p T > 0. Possible sign jumps of s T can occur for particular 
times where p T (t) = 0, and are decided by the normalization of the probability distribution. 

One may visualize this issue by a simple periodic time evolution for a two state system with 
t = 0, 1 [UJ. A periodic change between the probabilities po and p\ = 1 — pq can be described 
by the second order evolution equation d^po = 2w 2 (l — 2po), with solution po = cos 2 (u>i + a). 
In terms of the wave function this becomes a first order differential equation 

*(;)-(-SS)(s)- < 4 » 

The probability distribution is fixed for all t by the initial values po(^o) an d 5tPo(*o) f° r the 
second order evolution equation and equivalently by po(to), so(*o)*i(*o) f° r eq. @- Only the 
relative sign sosi matters, such that the needed "initial data" are the same for both descriptions. 
(It is actually possible to compute the signs s T (t) from "all time probabilities" [11 J.) For an 
expression of the relevant information contained in the wave function it is important that the 
latter is real. This contrasts with the discussion of classical mechanics in a Hilbert space by 
Koopman and von Neumann |12j . where phases of a complex wave function appear as new 
degrees of freedom beyond the probability distribution. 



2.2. Grassmann wave function 

A second useful formal tool for the description of our Ising type model is a map to an equivalent 
description in terms of Grassmann variables. This will make the close connection to fermions 
most transparent [10J. The Grassmann formulation relies on the isomorphism between states 
r and the basis elements g T of a real Grassmann algebra that can be constructed from the 
Grassmann variables ip 1 {x). Each basis element g T is a product of Grassmann variables 
g T = ip (x 1)^72 (5C2) • • • which is ordered in some convenient way. If a Grassmann element 
g T contains a given variable ip~,(x) we put the number n 7 (x) in the sequence r to 0, while we 
take n 7 (x) = 1 if the variable 4' 1 {x) does not appear in the product. 

An arbitrary element g of the Grassmann algebra can be expanded in terms of the basis 
elements 

9 = ^2qrgr- (5) 

T 

A time dependent wave function jg T (t)} can therefore be associated to a time dependent element 
of the Grassmann algebra g(t), provided the coefficients q T (t) are real and normalized according 
to ^r^r(^) = !• We will formulate the fundamental evolution law as an evolution law for the 
"Grassmann wave function" g(t). 

For this purpose we will formulate in the next section a Grassmann functional integral for a 
quantum field theory of massless Majorana fermions. It will involve -/V s L 3 /8 Grassmann variables 
ip^(t, x) for every discrete time point t. The Grassmann wave function g(t) and the fundamental 
evolution law can be extracted from this functional integral. This will establish a map between 
a quantum field theory for fermions and Ising type classical statistical ensembles with dynamics 
specified by an appropriate evolution law. From the point of view of the Ising type classical 
statistical ensemble the functional integral is a pure technical device for the specification of 
an evolution law for the classical wave function {q T (t)} and probability distribution {p T (t)}. 
Conceptually, we could start equivalently by specifying the matrix R Tp (t,t — e) in eq. ([3]). 
In practice, the functional integral formulation is much simpler since it provides for an easy 



realization of appropriate symmetries of the evolution law, as Lorentz symmetry in our particular 
model. 



3. Evolution law, action and functional integral 

3.1. Action and functional integral 

We start with the action S and Lagrangian L(t) for the regularized quantum field theory of 
massless Majorana spinors (N s = 4) 



Here we define 



tt-e 



t — tin 33 



B,(t + e,x) = l\ 



Y ~f6{{v})ip5{t + e,x k + v k A), 



(6) 



(7) 



with 



^(W) = o 1 ~ ^2( y k + w kI)T k - viv 2 v 3 I 



7# 



(8) 



Thus B~(t + e, x) involves a linear combination of Grassmann variables at lattice sites which are 
diagonal neighbors of x, corresponding to corners of a cube with basis length 2A, with x being 
at its center. The sums extend over j,k = 1 ... 3 and each corner corresponds to a particular 
combination of the three signs Vj = ±1. We define w k by W\ = —V2V^,W2 = t>it>3, W3 = —v\Vi- 
Eqs. ©-([H]) and the following imply a summation over repeated indices 7, 5 = 1 ... 4. 
We have also introduced the three real symmetric 4x4 matrices T k as 



0,1 
1,0 



To 



0,c 
-c,0 



1, 
0,-1 



0) 



where c = ZT2, r k are Pauli matrices and 1 stands for the unit 2x2 matrix. The product of these 
matrices yields the real antisymmetric 4x4 matrix 



c,0 
0,c 



TiT 2 T 3 , I 1 



-I. 



(10) 



The sum in eq. © extends over discrete time points t n , with f t = e = e , t n+ \ — t n = 
e,n £ Z, ti n < t n < tf. The time-continuum limit is taken as e — > for fixed ti n ,tf. Similarly, 
we sum in eq. © over points x of a lattice. For this purpose we consider for every given t a cubic 
lattice with lattice distance 2 A and J = 8 A 3 Yl x - ^ e ma y * a ^ e 6 = ^ suc h t na t the space-time 
lattice points belong to a hypercubic bcc lattice with distance 2A between nearest neighbors, 
which we call the "fundamental lattice". For even t = 2ne the space lattice points are even, 
x k = 2m' k A, with n, m' k £ 7L. This will be called the even sublattice. The odd sublattice consists 
of the odd time points t = (2ra + l)e for which the space points are also odd, x k = (2m' k + 1)A. 
The action ([6]) involves indeed only Grassmann variables for points on the fundamental lattice. 
For even t eq. ([6]) involves Grassmann variables tp-y(t, x) living on the even sublattice, while the 
combination B^(t + e,x) lives on the odd sublattice. For odd t the role of the sublattices is 
exchanged, now with tp on the odd and B on the even sublattice. This construction eliminates 
"lattice doublers" - a more detailed discussion of the lattice implementation can be found in ref. 

The action ([6]) is an element of a real Grassmann algebra - all coefficients multiplying 
the Grassmann variables ip~(t,x) are real. Within the Grassmann algebra the operation of 



transposition amounts to a total reordering of all Grassmann variables. The action ([6]) is 
antisymmetric under this operation, S T = — S. If we define formally the Minkowski action 
Sm = iS the latter is hermitean, Sm = S M since S* M = —Sm- 

The action ([6|) is Lorentz-symmetric in the continuum limit e = A, A — > 0. Indeed, performing 
the continuum limit the action ©-([E]) becomes [9] 

' {^dtip.y - ^{T^^sdkips}- (11) 

t,x 

It involves now four Grassmann functions ^~ f {t,x),^ = 1 . . . 4, x = (xi,X2,xs). The integral 
extends over three dimensional space and time, with dt = d/dt and dk = d/dxk- The Lorentz 
invariance of the action ([lip is most easily established by employing the real matrices 

y°={-° ; n „).7* = -7% (12) 

such that 

S = - [ tfW > $ = V'V, (13) 

Jt,x 

where ji = (0,k) and do = dt- The real 4x4 Dirac matrices 7^ obey the Clifford algebra, 
{7^,7^} = 2r]^ u , with signature of the metric given by r]^ u = diag(—l, 1, 1, 1). 
The functional integral is defined by the partition function 

Z = j V^g f [i;(t f )}e- S g in [^(t in )] , J^ = Y[J (#4^, x) . . . #i(t, x)). (14) 

The boundary terms gi n and gj only depend on the Grassmann variables ip(ti n ) and ift(tf), 
respectively. The boundary terms g~f and gj n are related to each other [9], such that the functional 
integral (|14p is fully specified by the choice of gi n . In turn, gj n is given by the initial values for 
the classical wave function, cf. eq. ([5]). 

3. 2. Grassmann wave function from functional integral 

In this section we compute for the functional integral (j!4[) a Grassmann wave function g(t), which 
is an element of the Grassmann algebra constructed from the Grassmann variables ip^(t,x) at 
given t. The functional integral (j!4|) involves variables for arbitrary time points t n . In order to 
construct a wave function g(t) which only refers to a particular time t we have to "integrate out" 
the information referring to other time points t' 7^ t \10\ 113]. This can be done by decomposing 
the action ^ 

S = S < + S > , S < = J2 L (t'), S > = J2 L (t')- (15) 

t'<t t>>t 

The wave function g(t) obtains now by integrating out all Grassmann variables for t' < t 

g(t) = J Vi>(t' < t)e~ s <g m . (16) 
As a consequence of its definition (|16p the Grassmann wave function obeys the time evolution 
g(t + €)=/" Vi,(t)e- L{t) 9(t) = [ Vi/}(t) exp { - ^ ^ ^(t, x)B y (t + e, x)}g{t). (17) 

J J x 7 



This determines {q T (t + e)} in terms of {qy(i)}. Thus the action ©-([8]) specifies the dynamics 
how the probability distribution {p T (t)} evolves in time. The particular dynamics of a given 
model is determined by the form of i? 7 in eq. (JSJ). In a similar way we can relate g(t + 2e) to 
git). Performing the necessary Grassmann integrals over ip(t) and t/j(t + e) one obtains [9j for 
e — > the differential evolution equation 

digit) = j- e [git + 2e) - g(t)] = Kg{t) = ^ Qr(t)>Cg T = ^ q T {t)g p K pr = ^ d t qAt)g T , (18) 

T T,p T 

with Grassmann evolution generator K, 

K = T,d^^d k ux). (19) 

From eq. (I18p we can directly extract the evolution equation for the classical wave function ([3]) 
in a continuous version, 

d t qr{t) = ^K Tp q p it) , K pT = I V^~g p Kg T . (20) 
p J 

We identify this evolution equation with the generalized Schrodinger equation for a quantum 
wave function for the special case of a real wave function and purely imaginary and hermitean 
Hamiltonian H = ihK. It specifies the evolution of the probabilities for the Ising-type classical 
statistical ensemble. 



4. Majorana fermions 

4-1- Conserved quantities 

For a physical interpretation of the evolution equation (fTgj) , (fH?|) it is useful to look for conserved 
quantities. They correspond to time independent expectation values of classical statistical 
observables. We first represent the classical observables in the Grassmann formalism. Consider 
the observable measuring the occupation number iV 7 (x) of the bit 7 located at x. The spectrum 
of possible outcomes of measurements consists of values 1 or 0, depending if a given state 
r = [n 7 (x)] has this particular bit occupied or empty. For the Grassmann basis element g T 
associated to r one finds N-yix) = if g T contains a factor tp-y(x), and N y (x) = 1 otherwise. We 
can associate to this observable a Grassmann operator A/" 7 (x) obeying (no summation here) 

d 

AA 7 (x) 5t = (N y (x)) T g T , N y (x) = — j—^ ifij(x). (21) 

Two occupation number operators A/" 7l (xi) and A/" 72 (x2) commute. 

In general, we may associate to each classical observable A a diagonal quantum operator A 
acting on the wave function, defined by iAq) T = A T q T with A a diagonal operator A rp = A T 5 Tp . 
In the Grassmann formulation one uses the associated Grassmann operator A obeying 

A 9r = A rg7 , { A ) = Jv r9 A 9 . (22) 

Here g is conjugate to g, i.e. for g = J2 T Qt9t one has g = ^ T q T g T . The conjugate basis elements 
of the Grassmann algebra g T are defined [10] by the relation 



3r5r = |0)=nil^) > =1, 

x 7 J 



(23) 



(no sum over r) and the requirement that no variable ip^{x) appears both in g T and g T . 

Using dtq = Kq (j20|) and dtg = ICg (fl~8|) the time evolution of the expectation value is 
expressed in terms of the commutators 

d t (A) = (q[A,K]q) = f Vi;~g[A,IC}g. (24) 

Conserved quantities are represented by Grassmann operators that commute with /C, [A, K] = 0. 



4-2. Particle states 

Our system admits a conserved particle number, corresponding to the sum over all occupation 
numbers iV 7 (x) in the lattice, or to the Grassmann operator AT, 



d_ 



M = I 7777777^) > W, tC] = 0. (25) 



We can decompose an arbitrary Grassmann element into eigenstates of J\f, g = 
Mg m = mg m . For an eigenstate g m the probabilities in the Ising-type classical 
statistical ensemble differ from zero only if the total number of "occupied bits" equals precisely 
m. The time evolution does not mix sectors with different particle number m, such that the 
coefficients A m are time independent, dtg = T, m A m dtg m , Qt9m = K-9m- We can restrict our 
discussion to eigenstates of M. 

A static vacuum state go with a fixed particle number mo obeys 

Kgo = , Mgo = m g , J V^g g = 1. (26) 

An example for a possible vacuum state is the totally empty state go = \0) (eq. (|2"5]1). Another 
example is the totally occupied state with go = 1, mo = B = N s L 3 /8. There are many possible 
vacuum states, for example a half filled state with mo = B/2. 

We next define creation and and annihilation operators a\(x), a~(x) as 

d f 
a\{x)g = — j-^ g , a^(x)g = ip y (x)g , {a\(x), a e (y)} = 6 je 5(x-y) , N = J aj f (x)a 7 (x). (27) 

Acting with the creation operator on the vacuum produces one-particle states 

9i(t) = / q 1 (t,x)a\(x)g = Gi9o , (AA - m )gi = g\. (28) 

J X 

The time evolution of the one-particle wave function g 7 is given by 









IX 


JC — 

' dtp 


9o = 

J X 



j>4 )s °- (29) 



dtgi = fcgi 

Inserting the specific form of the evolution operator p9p yields the Dirac equations for a real 
wave function |1U] 

7"0 M g = O. (30) 

We emphasize that the Dirac equation follows for arbitrary static states go which obey /C<?o = 0. 
Since there is no mass term and the Dirac-matrices as well as the wave function are real, eq. 
(|30p describes the time evolution of the wave function for a single massless Majorana fermion. 



States with of n fermions can be constructed by applying n creation operators a\(x) on the 
vacuum. For example, the two-fermion state obeys 

92(t) = -j= I q J£ (t,x,y)a\(x)at(y)g . (31) 

V ^ J x,y 

Since the creation operators anticommute the two-particle wave function is antisymmetric, as 
appropriate for fermions. From there we recover the standard interference rules for identical 
fermions. The linearity of the evolution equation in q is crucial in this respect. It guarantees 
the superposition principle for the wave function, and therefore interference effects for the 
probabilities. 

5. Dirac fermions in electromagnetic fields 

Dirac spinors can be composed of two different Majorana spinors with equal mass (N s = 8,B = 
L 3 ). We extend our model by introducing in the action a mass term ~ m and a coupling ~ e to 
an external electromagnetic potential A^(x), (with suppressed spinor indices) 



{Mdt ~ T k d k + m 7 °/)Vi + Mdt - T k d k + my°I)i> 2 } 

{V>i(A) - A k T k )ip 2 - ~ A k T k )ih}. (32) 



t.x 



We introduce a complex structure by defining the four-component complex Grassmann 
variables ipo = ipi +iip2 , V'z? = V'l — ^2- In terms oiipD the action (j32[) takes the familiar form 

$D('fd f j, - mI)ip D -ie tjj D ^A^ip D , ip D = ^7°- (33) 

t,x Jt,x 

5. 1 . Evolution equation for Dirac fermions 

For the action ([32]) the Grassmann evolution equation dtg = K,g involves K. = /Co + JC m + K-a 



)C + tC m = / 77TT -r(Tkd k -m 7 °I)Mx), (34) 

Jx a=l,2 a ^ a y X > 



JCa 



-(A (x) - A k (x)T k )^ 2 {x) - ——(A (x) - A k {x)T k )Mx) 



dipi(x) difi2{x) 



It describes the dynamics for an arbitrary number of charged relativistic fermions (and their 
antiparticles) in external electromagnetic fields. 

The evolution equation for the classical wave function (|20p . and therefore for the probability 
distribution, obtains again from eq. (|18|) . In summary, we have formulated a time evolution 
equation for a classical wave function q T (t) which describes an arbitrary number of charged 
electrons in external electric and magnetic fields. Its restriction to one-particle states yields the 
relativistic Dirac equation. As usual, a non-relativistic approximation will yield the familiar 
Schrodinger equation for a particle in a potential or moving in external magnetic fields. 



5.2. Classical observables 

Classical observables can be constructed from linear combinations of products of occupation 
numbers. We combine the index a = 1,2 for the two Majorana spinors ip a with the index 7 



into a common index e = (7, a),e = 1 . . . 8. The Grassmann operators corresponding to the 
occupation numbers N e (x) obey eq. (|21|) . A local particle number can be defined as 



N{x) = J2K(x) , [M(x),K] = (T k ) m d k M ari , M ev (x) = ^-j-^ v ( x ). (35) 

According to eq. (|24|) the time evolution of the classical mean local particle number obeys 

d t (N(x)} = (T k ) va d k (M ari ). (36) 

This involves the expectation values of "off-diagonal" operators .M Qr) . Products of operators 
M QT j can be defined by the non-commutative quantum product and we see here a first appearance 
of non-commutative structures. The total particle number J\f = J N{x) , [J\f, K] = 0, is 
conserved. 

6. Quantum mechanics for particle in a potential 

6.1. Dirac equation 

We next concentrate on one-particle states, described by Grassmann elements 

9i{t) = L \ qi ^ x) d^ + q2 ^ x) d^) 9 °- (37) 

Here go is some arbitrary static "vacuum state" . For ICgo = we find for the one-particle wave 
function the Dirac equation 

7 /i (<^ + ieA^ipn = im'yipD , ip D = qi + iq 2 . (38) 

This equation holds for an arbitrary representation of the Dirac matrices 7^. The matrix 7 
(often called 7 s ) obeys 7 = — i7°7 1 7 2 7 3 . In our conventions [T3| the Dirac equation obeys 
parity, charge conjugation and time reversal symmetry. 

The derivation of the Dirac equation for a one-particle state (|37p has only used the general 
evolution equation for Grassmann elements g(t) and the condition /C<?o = 0. It therefore describes 
the dynamics for a very extended family of classical probability distributions or classical wave 
functions. Indeed, it is sufficient that go is an arbitrary static state (not necessarily a priori with 
a fixed particle number). This reflects the physical property that isolated one-particle states 
can occur under a wide variety of circumstances, and that for sufficient isolation the properties 
of the environment do not matter for the dynamics of the isolated particle. In this context we 
emphasize, however, that our model only describes external electromagnetic fields while we do 
not account for the fields generated by the particles that may be present in the environment. In 
this sense our setting describes "real physics" only in a situation where the electromagnetic fields 
generated by all present particles can be approximated by a "mean field" that is independent 
of individual particle positions. This is precisely the setting of standard one-particle quantum 
mechanics. 

6.2. Particle observables 

Let us consider the expectation value of the local occupation number M{x) defined by eq. ([3"5]) , 
For go = |0) its expectation value obeys the simple relation 

(at(x)> = J2^ x ) > / E M*)) 2 = / rf>(*)M*) = i- (39) 

e J X £ J X 



The normalization (j39|) is preserved by the unitary time evolution of the one-particle wave 
function. Indeed, we can write the Dirac equation (j38|) as a Schrodinger-type equation 
ihdt^pD = Hipjj, with hermitean Hamiltonian H = ihT^dk + hmj ^ + fte(Ao — T^Afc). 

For a pure one-particle state the expectation value (N(x)) can be interpreted in a natural way 
as the probability density to find the particle at the position x. This is precisely the standard 
interpretation of i/)' D (x)t/j d(%) in one-particle quantum mechanics, 

w{x) = (N(x)) = ^ D (x)ip D {x) , / w{x) = 1. (40) 

J x 

Since N(x) is a classical observable, we can define the position of the particle as a classical 
observable, and similarly for functions of X, 



X = / xN(x) , f(X) = / f(x)N(x). (41) 

J X J X 

The expectation value in classical statistics coincides with the standard quantum mechanics rule 
(X) = ( f xN(x)) = [ xw(x) = [ ^ D (x)xip D (x) , (f(X)) = [ ^ D (x)f(x)cp D (x). (42) 



In particular, one obtains the same formula for the dispersion (X^X^) — {Xk){Xk) as in quantum 
mechanics. We conclude that measurements of the position of a particle, or more generally the 
distribution of positions in an ensemble of one-particle states, can be described equivalently in 
an Ising type a classical statistical ensemble with classical observables, or in quantum mechanics. 
The complete time evolution of the distribution of positions is identical in both descriptions. 
This covers, in particular, the characteristic quantum interference in a double slit experiment. 
The time evolution (|20p . (|38p for the classical wave function and associated classical probability 
distribution {p T } produces for one-particle states exactly the quantum mechanical interference 
pattern. 



6.3. Schrddinger equation 

Standard quantum mechanics for an electron in a potential is recovered from the non-relativistic 
approximation to the Dirac equation. This is well known, and we sketch here for completeness 
only the case = 0. The nonrelativistic approximation becomes valid if cAq and iT^dk 
are small compared to m. Since (7°7) 2 = 1, it is convenient to choose a new basis where 
7°7 = diag(l, 1, —1, —1). With V(x) = heAo(x) and M = Km the Hamiltonian takes the form 

*=(4L : _M l ) +l ' w ' pi= -* M " (43) 

with matrices o~\ = — 1 , 02 = IT2 , 03 = 1T3. The appearance of h only fixes the units for pk,M 
and V(x). 

We next decompose (pr> into two-component wave functions, (p^ = {\ T ,p T ). For the non- 
relativistic electron we consider the approximate solution p = Ax, where A = a\pk/{2M) 
is determined by requiring in leading order dtp = AdtX- Insertion into eq. (|38p yields for 
ip = exp(iMt/H)x the standard Schrddinger equation for a particle in a potential V, 

= {^w + v ) ( 44 ) 

All quantum mechanical phenomena extracted from solutions of the Schrddinger equation 
are described by our time evolution equation for a classical statistical ensemble of Ising-spins on 



a lattice. For a potential realizing a double-slit situation the standard interference pattern will 
appear behind the slits. This holds provided that the initial state at some time to corresponds 
to a one-particle state describing a particle moving towards the slits. Interference is realized by 
classical probabilities. 

For the choice go = |0) the evolution law for the classical probabilities associated to a one 
particle state is rather simple. The only non- vanishing probabilities p e (x) = q^(x) occur if 
precisely one lattice site is occupied with a given species e. The evolution of q e (t,x) obeys eq. 
(|38p . We recall, however, that the same physical situation arises for arbitrary static go- The 
time evolution of probabilities can then get rather complicated in the Ising-spin picture whereas 
it remains simple in the fermionic picture. 



6.4- Particle-wave duality 

The discreteness of measurement values in quantum mechanics can be traced back to the discrete 
occupation numbers of the Ising-type model. In quantum mechanics we may define an "interval 
observable" J-ji by a function 

J ^ = { I otherwise ) ' J * = ^ ^ 

such that its spectrum consists of the discrete values and 1. According to the rules of quantum 
mechanics, the possible outcomes of a measurement of J-ji are or 1. The interpretation in 
quantum mechanics is simple: either the particle is within the region (interval) Jn, in which 
case the measurement value Jn = 1 will be found, or it is outside this region, and Jn = will 
be found. Particles are discrete objects - they are either inside or outside an interval. 
In our classical statistical Ising-type setting Jn is a classical observable, given by 

J n = f N(x) = Y^N L (x). (46) 
Jn n 

The sum Yin extends over all lattice points within the region 1Z. The classical observable 
Nl(x) corresponds to a normalization of occupation numbers for the discrete lattice where 
{N L ,e) T = 0, 1. The sum over species Nl(x) = Yl NL, e (x), cf. eq. (f35|) . can therefore take in any 

e 

classical state r only the discrete values (Nl(x)) t = (0, 1 . . . , 8), according to the total number 
of eight species. In consequence, for any state r of the classical statistical ensemble, (Jn)r is 
a positive integer or zero. According to the standard rule of classical statistics these integers 
describe the possible outcomes of measurements. 

For a one particle state the total particle number equals one, 

1 = J N(x) = Y,Nl(x), (47) 



V 



where Yl extends now over all lattice points in the total volume. Since 1Z must be contained 

V 

in V the maximal allowed value for Jn in a one-particle state is one, such that {J-r) t = 0, 1 
are the only possible values of the classical observable for such a state. This is precisely the 
quantum mechanical rule. No new postulate is necessary for this measurement in quantum 
mechanics - the quantum rule is inferred from the standard rule of classical statistics. Of course, 
the expectation value of J-ji for a one particle state is the same in the quantum mechanical and 
the classical statistical description 

(Jn) = / (p f D (x)J n (x)(p D (x) = / ip\j(x)<pr>(x). (48) 
Jx Jn 



While our model connects the discrete particle aspects in quantum mechanics directly to 
the discrete classical Ising-spins, the continuous wave aspects also arise in a natural way. 
The quantum wave function is continuous because probability distributions and associated 
classical wave functions are continuous (at least piecewise). As we have mentioned already, 
the characteristic interference effects for waves and the superposition arise from the linearity of 
the fundamental evolution equation in the classical wave function. 

7. Conclusions 

We have derived the Schrodinger equation for a quantum particle in a potential from a classical 
statistical ensemble for Ising-spins. The dynamics of the classical statistical system has to be 
specified by an appropriate evolution equation for the probability distribution. For this purpose 
we employ the classical wave function which is defined as the positive or negative root of the 
probability distribution. The proposed evolution equation is a linear differential equation for 
the classical wave function. The wave function at time t obtains from the wave function at t' by 
a rotation - this guarantees the preservation of the normalization of the probability distribution. 

We have exploited a map between the classical wave function and a Grassmann wave function 
which is an element of a real Grassmann algebra. In turn, the time evolution of the Grassmann 
wave function can be associated to a Grassmann functional integral. This allows us to formulate 
the evolution equation for the classical wave function in terms of the action of a functional 
integral. The symmetries of the model, as Lorentz symmetry and electromagnetic gauge 
symmetry for our model of Dirac spinors in an external electromagnetic field, can be easily 
implemented in this way. Since the classical wave function is real we have to formulate the 
model in terms of a real Grassmann algebra. 

Our model is regularized on a lattice of space points. On the one hand, this guarantees that 
mathematical expressions are well defined for a finite number of lattice points, with continuum 
limit of an infinite number of points taken at the end. On the other hand, the concept of classical 
Ising-spins or associated occupation numbers at every point x is well defined. The discreteness 
of the particle aspects of quantum mechanics can be traced back to the discrete occupation 
numbers that can only take the values zero or one. We also have used discrete time steps, 
and we have formulated the lattice action such that the Grassmann wave function g{t + 2e) 
obtains from g{t) by a rotation. The time evolution is unitary not only in the limit e — > 0, but 
also for finite e. A unitary time evolution for infinitesimal time steps is easily achieved by any 
antisymmetric evolution generator K Tp for the wave function, q T (t), cf. eq. (|20p. Generating 
a unitary evolution also for finite smallest time steps imposes restrictions on the form of the 
lattice action. This, together with the requirement of a real Grassmann action, requires some 
care for the construction of the model and explains the specific form of the action ©-([H]) in 
comparison with other possible lattice actions. 

The proposed evolution equation for the classical statistical ensemble of Ising-spins does not 
only lead to the Schrodinger equation for non-relativistic one-particle states. It entails the full 
dynamical equations for a quantum field theory of Dirac fermions in an external electromagnetic 
field. The dynamics of states with an arbitrary number of fermions, including the characteristic 
interference patterns for indistinguishable fermions in quantum mechanics, is correctly described. 

At this point it seems worthwhile to ask some questions about the origin of characteristic 
features of quantum mechanics in our classical statistical setting. Particle- wave duality is realized 
by the discreteness of Ising-spins on one side, and the continuous probability distribution or 
classical wave function on the other side. Interference arises from the formulation of the basic 
evolution law in terms of the classical wave function. While the classical wave function is real, 
it can nevertheless take positive and negative values which can add to zero locally if two wave 
functions are added. The superposition principle or linearity of the quantum evolution finds a 
direct origin in the formulation of a dynamical law for the classical statistical ensemble that is 



linear in the classical wave function. The characteristic physics of phases in quantum mechanics 
is connected to the presence of a complex structure within the real Grassmann algebra. Planck's 
constant h appears purely as a conversion factor of units. The uncertainty relations can be 
obtained directly from the possible solutions of the Schrodinger equation. 

Finally, one of the most characteristic elements of the mathematical formulation of 
quantum mechanics is the presence of a non-commutative product for operators and associated 
observables. These structures are obviously present in our formulation and clearly very useful 
for a discussion of solutions of the Schrodinger equation or Dirac equation. The essence of the 
emergence of non-commutative structures is the coarse graining of information [TS] , [3] , [2] . 

This issue has not been addressed in the present note. It seems reasonable to expect, however, 
that a system that is governed by the Dirac equation on microphysical scales - say lattice 
distances shorter than the Planck length - will also show similar properties at "macroscopic 
scales" associated to coarse graining. (Such macroscopic scales can still be much smaller than 
all characteristic scales of atom physics or elementary particle physics.) The basic reason is that 
the form of the Dirac equation for one-particle states is essentially fixed by the symmetries. It 
will not be altered if the coarse graining respects the symmetries. While it remains an interesting 
task to perform this coarse graining explicitly, the main message of this note is already very 
clear at the present stage: Quantum field theory and quantum particles can be obtained from a 
classical statistical ensemble. 
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